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I . INTRODUCTION 


IC hab been recently recognized that aubmllllmeter waves are most desirable 
for remote sensing in tha atmosphere. The interest in the generation of coherent 
submilllmeter waves has made the study of superradiance under realistic conditions 
an urgent need [1]. Ihe central problem in this study is the Interaction of multi- 
mode radiation with extended resonant media. Using the extremely powerful, but not 
very well known, Dlcke-Schwendimann representation [2,3], we should be able to put 
the coupled superradiance master equations, recently derived by Picard and Willis 
[A], in an explicit form. With this as a starting point, we can very easily study 
the effects of superradiance on pulse propagation in resonant media. 

In thb period of thirteen months. Wo have had very significant achievements 
toward this goal. However, the whole project is too big to be finished in one year. 

We will describe the major results of our studios. The description will be brief 
since the details of these results has been published or to be be published in open 
literature. 

II. MATRIX ELEMENTS OF COLLECTIVE OPERATORS 

One of the chief difficulties encountered in the Dlcke-Schwendlraunn approach 
is the lack of explicit expressions for the matrix elements of collective operators. 

We hav'* developed a simple and systematic diagrammatic teclmlque to clear up this 
difficulty. The detailed description of this technique has been published [5J. 

III. TRANSITION FROM INCOHERENT TO COHERENT RADIATION 

In the collective emission of radiation, the cooperation number r of the atomic 
system, as defined by Dicke [2], can remain the same; out it can also change. The 
i.atter is often ignored. We have found that the emission through r-conservlng process 
is coherent and that through r-nonconserving process is Incoherent. The uransitlon 
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from Incoherent to coherent radiation can be beat described by conaiderlnR the 
number of photons emitted, instead of the time, as the Independent variable. 

The former is a discrete variable while the latter is a continuous variable. 
Hence, we have a difference equation instead of a differential equation. The 
result of the numerical solution of this difference equation has been described 
in detail in a paper to be published [6]. 


IV. RATE EQUATION FOR SUPERRADIANCE 

The difference equation described in the p* . ms section is not suitable 
for the study of the intensity of superradiant pulse. For this purpose, we have 
obtained a rate equation as follows: 


dP(r.m) _ (r-fpH-2) (r-Hn-fl) p(ri-l.m-fl) + (r+m+l) (r-m)vP(r ,ra+l) 
dt 2rf2 

(rfm) (r-Hn-l)p (r.m) - (r+m) (r-m+l)YP(r,m) , 

2r 

where r is the cooperation number, ro is one-half of the difference of the numbers 
of atoms in the excited state and the ground state, P(r,m) is the probability of 
the atomic system in the Dlcke state |r,m>, and y is the radiation rate of a 
single atom. 

It is convenient to make the following change of variables: 
k ■ r-Hn, t, ■ r-m; 

and to define the following generating function: 

N N 

f(x,y,t) • 'I I P(k,l)x*V*^* 
k-1 i-1 

Then we can replace the rate equation for P(r,m), a function of t, by a partial 
differential equation for f(x,y,t) as follows: 
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The solution to this last equation is yet to be found. 
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IV. PUBLICATIONS 

Two papf>m have been prepaired based on the results of this raaearch project: 

1. ''DlnKfsinnaClc Technique for Calculating Matrix Elenencs of Collective 
Operators In Supcrradlance" 

Physical Review A| Vol. 12, No. 2, pp. 575-386, August 1975. 

2. "Difference Equation for Superradiance" 

To be published.' 
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DiaKomnuitic technique for calculatinK matrix elementx of collective operators in 

xuperradiance* 

China T%ung |4W 

OrpurlmrHt of Phy\n\ and Uaihrmolit x. 4lohomo 44 W I'niKrmiy. .Marmot. Alohtmo tiftj 

(KKcivfd 10 April l«7)| 

Ailopling ihc MvcaJIcd *‘|m<al<iti«.al vontlriKlMm.'* one can rtprcM the ri(.nMalc» of collcclivr oprralm 
corrnpondini In ■ cpcciflct* motlr for an .V-alom cyclnn in Irrmt of ihow for an I \ ly-altan cyalcm I'tinf 
ihcar l>H'k« Main at hatn and uting Ihc Wigncr-I ckail IhriHrm. a mama ckmcnl of a coliccnvr ojirriliH of 
an arhiirary mode can he wnllen aa the produci of an m drpendenl faclor and an m independeni reduced 
malm elrmenl iRMKi. A wl of rccuranin rormulct for the KMI la obtained A fraphnal repreacnlalMin of 
Ihc KMI on Ihc braiK'hing diagram for hmary irreducible reprcamlalHaia id permutalum groupa m Ihen 
iniroduced Ihia givn a Mmpie and ayalemalic way of cakulaling Ihc KMI’ Thia method it eapecially taeful 
• hen ihc ctaiperalion number r la cloae lo .V/2, ahere almoal caaci aaymplofic eapreaaiona can he oblai;«d 
eaail) The mull ahoaa eaplKilly the geometry depcndeiwe of aupcrradiaiHC and ihc rdaliae importaiKc u.' r- 
conterving and r ntmconterving proccaan Ihia cicara up the chief difficulty nwouniered in ihe Dnlie 
Schaendimann approach lo Ihe problem of V lao leaH aloma. apread over large reguma. inlcracting aiih a 
mullimodc radialHMi field 


I IMKOIIU I ION 

In 1954 Dirkp' diacung(<d Ihe spuntaneoun emiii* 

■ lun of radiation from an excited ayatem of N 
identical two-level atoma. Cy conaiderint; the 
entire collection of atoma as a ainKle quantum- 
mechanical ayatem, he found tliat under certain 
conditions the Individual atoms cixiperate to emit 
radiation at a rate proportional to ,V', which Is 
much greater than their Incoherent emission rate. 
This phenomenon Is called super rad lance. 

In Ihe vast literature on superradlance, ' n.jst 
of the Investigations have been limited to cases 
in which the emitting atoms either are confined to 
a region smaller than the wavelength of emitted 
light or are able to couple with only one radiation 
mcxie. This situation permits important mathe- 
matical simplifications. However, current at- 
tention’'* in this field is being focused on the more 
difficult, but more realistic, problem of A two- 
level atoms, distributed over a space of dimension 
much greater than the radiation wavelength, inter- 
acting with all the mtxles of radiation field. Among 
the works devoted to this problem, the one by 
Schwendimann* elaborates a very elegant approach 
of the problem initiated by Dicke.* One of the 
chief difficulties encountered in this approach is 
the lack of explicit expressions for the matrix 
elements of collective operators. In this paper, 
we will develop a simple and systematic diagram- 
matic technique to clear up this difficulty. 

II VtODIt H\Mlt10MY> 

We consider a system of ,V identical two-level 
atoms distributed over a space of dimension much 
larger than the wavelengths of the multimode 


radiation with which it interacts. In the usual 
electric dipole and rotating-wave approximations, 
the system ran be described by the following 
Dicke Hamiltonian; 

II /i^ wf rtfuf ^ If 


where k is the wave vector, oj and ut are the 
corresponding photon creation and annihilation 
operators and obey the usual bision commutation 
rules, and is the corresponding frequency, 
w, IS the transition frequency between the two 
levels of each atom. H':, and H, are collective 
atomic operators defined in terms of the single- 
atom flip operators H',, H,, and /<„ , which •'S^v 
angular-momentum commutation rules as ' 

“ f I 

where is the position of the ith atom, and 
>ff ■(•*'7 2i-or*)‘'p„ are the coupling constants, 
where p„ is the transition dipole moment, and 
V* is the volume of Ihe cavity in which the radia- 
tion field IS quantised. 

As one sees from Flq. (2), th** Hi operators de- 
pend on the wave vector k and jn Ihe positions 
}{,'s of the atoms. Their commutation relations 
are 

(«>./<:. I . (3) 

I I 

Here we see that the commutator of two opera- 
tors corresponding Indifferent modes given an 
operator corresponding to a third mode. There- 
fore, in contrast to Ihe case of samples smaller 
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than a wavelpngth in which Uir ' (artcira dis> 
appear, the behavior ot the atoms cannot be des> 
cribed by a alniple set of an|{ular>momentum op- 
erators. However, It la possible to factor out 
the C dependence In the operators and confine It 
to r-number functions. This la the key point of 
the Dlcke-Schwendlmann approach, to be dis- 
cussed In the folluwlmt. 


K Uskr VhsmdMnMin ippriiat h 

Consider the set of operators |/f£|, as 

defined by fCqs. (2) with a particular wave vector 
C, which we will call the principal mtxle. It can 
be easily shown that they obey the anKular-momen- 
tum commut|tlm relations. So we can chose 
I ♦), and to be the excited and ((round 

states of the ith atom. Then we can construct 
eixenstates |r,Mt; a) of and A, sucl. that 


J r. a) • r(r ♦ 1)1 r, *«; a>. 


(4) 


where ■ J and a ■ I, 2, . . . 

IS a degeneracy iniex with 


,* ^’U2> ♦!) 

' ' (,V 2*r ♦l)l(^/a-» )l• 


(5) 


The eigenstates |r, *H;ai) will be called Dicke 
states, r will be called the ccMiperatlon nunibe;*, 
and in Indicates the population Inversion or t<>tal 
energy of the atomic system. The trouble Is that 
the eigenstates of corresponding to the prin- 
cipal nuxle are not eigenstates of h\ corn*spond- 
ing to ot.ier modes. So we might think that we 
need to construct a set of Dicke states for each 
ptissible mode, which would be an imposing task. 
Fortunately, this is unnecessary. 

If we define’ 


«i:-i (6> 

. I 

(or an arbitrary mode k and consider the set of 
operators {«:, H*. A,;}, they do not obey 


anguUr-.nommtum commutattun rules. H'lwever, 
as first pointed out by Dicke,' their commutation 
relations with the established anf ular-momentum 
operators and H, are of the following form; 

J’l. -.3. (7) 

where we should Identify ff*. /C. and H,; with 
T*., r'., and 7^, respectively. This proves* that 
)/('*, if', ff,;} IS an irreducible tensor operator 
set of rank 1, or simply s vector operstor. As 
a consequence, we esn discuss them in terms of 
the eigenstates \ r,m\a) of /fr, operators. Fur- 
thermore, the Wigner-Eckarl theorem can be 
used to establish the selection rules for the ma- 
trix elements (r' , w'; r, w; o>; namely, 

r-r'stlorO and m~m' ■ t I ; and the nonvan- 
ishing matrix elements can be factoriced as fol- 
lows: 

{r.mt l;o'l/f!:l»'. 

■ |(r <»N)(r tm* I)|'^(r; o'|/fr|r, o), 

( r ♦ t , »« 1 1; »*’ I /fH ♦'t »» : **) 

• » |(r t »»» t)(r t m»2)|''^(r » I; o'|7fh 


^ r - I, I; u'|/(V| r, M«; o) 

(H) 

where the reduced matrix elements (RMK) are 
independent of m and are of the following forms; 

(r; I A* I r; o) . ^ | c •'< * 

I 1 


(9) 


(r , I; lAfh »•: «») • t c;*r, c -'f 


■ ^ 


• on 

99 l| 


(k-k,). 


H Triiti%fomtr«l Hamillotiuii 
We now expand H\ In termn of Dicke BtateH 


r M a 9 m a 

• (C - k.MK ‘c*’ * (J - k.)(/”c., ♦ i. 


UO) 


where we have used Lqn. (8) and (9) and have defined 
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(«'C ^1 (*■♦.«)(* t m * I r. **• i I; 1 *')^ » , w; II 1 , 

m 

l/'*C" * *^|(»‘ » '»»♦ l)(r » »w ♦2>|‘^| » ♦ l.Mii I; o'H I", »m; <»l , (lU 

l)|‘*|r- l.m* 1; a'>< r. w; <i| . 


UilnR thrse rraulUi, the HanilRonun In ICq. (1) ran be rewritten aa 


iifaf 

* 

♦;r-'(C -5,)(/rC"’ ♦ h . <» 2 ) 

♦ H e.) ♦ (C-k.M/*T.*:, ♦ H e.)] . 

a 

where dearrlbea proreaaea In which the cooperation number r la coiiaerved, and deacribea pro* 


ceaaea In which r la not conaerved. 

Ill IVMl AtHINOI HIIHI III M<UKIX IIIMINIS 

In the Dlcke-Schwendimann approach that we 
have elaborated In Sec. II. we notice the impor- 
tant role played by the reduced matrix elementa 
which we have written aa ♦^""*(k-k,L a r-number 
function of the wave vectora. All the geometry and 
mode dependence of our problem are determined 
by theae functiona. They are independent of time, 
and hence, deacribe the atatic effecta of the aya* 
tern. However, theae functiona appear aa ci ef* 
firienta of dynamical uperatora In the expa'.alon 
of the Hamiltonian. Therefore, the uaefulneaa of 
thia approach dt*p«-iKla on whether we can obtain 
an explicit expreaalon for t*ach of them. 

\ '(•rnriliifH-al iiimlnk Imni” 

From the group-theor»*tical point of view, the 
aet of the ao*called avmmetrized atatea' |r,m:ut.. 


o a I, 2 j! forma a baaia for an irreducible 

reprt'aentation of the permutation group cor* 
reapiHidIng to a certain inary partition of .V. 
There la a one- to*one correapondence between a 
binary Irreducible repreaentation and the cooper* 
ation number r. Recauae the matrix repr'aenta* 
tion of a group la fixed only within a aimilaritv 
tranaformation. the explicit exprf*aaiona for RMK 
depend on how the Dicke atatea are conatructed. 

So we niuat choae a definite way to cunatruct 
them A very natural and convenient acheme la 
the Bo*called ‘‘genealogical conatruction.”* We 
will uae It here. 

In thia ronatrucllon, the eigenalalea of ,V alonia 
with elgenvaluea (•',»<) are conatructed from thoae 
of .V-1 atoma with eigenvaluoa (r - i) and 

(r-i, m* i), or (r J) and (r ♦ J, w* i). by 

adding the “apln“ functiona for the >Vth atom aa 
followa; 


|.V, r.w, o> •((»■ ♦»«) 2» |‘^| A - 1.*' - j. w- 1; «)l ♦;» ♦!(»■ -'«) 2* l.r - o, | e*‘‘'i‘'». 

for a ■ 1 , 2 f f 

|.V,r.m;a).|(»--»«*l) (2r ♦ 2)|‘^| .V - 1. r - i. m- i ; *1)1 ♦ )» 

-l(r emeU (2r ♦2)|‘*|.V -1. re J.w *1:3)1 

for i» «f *3 and 3" 1. 2 


and we have the relation 

which can be eaailv verified with the uae of Eq. 
(5). 

Starting from the elgenatatea for one atom, we 
can conatruct the eigenatatea for S atoma bv re* 
pealed application of Eqa. (13). The Dicke atatea 


olMained in thia way form an orthonormal aet 
automatically In thia conatruction, an eigenatate 
la apecified, in addition to r and m, by a aeriea 

of ‘partial*cooperation numbera" r,.r, r, 

which will replace 'he aingle index u. So we will 
adopt a new notation 

I V, r . m; <») a 1 r. »«; r,r, • • • r* > ( H) 
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for the piKFiMUtM. is alwaya i; r, ran br 
plther 0 or 1; tf r, la 0, r, can only bp i, but if 
r, la I, r, ran bp pithpr i or 1; and ao on. Anv 
way. wp muat havp 

r, 0 and mt (15) 

H Rn ur\KNi liirmuU Uh mlutni mjin« rlrmrnlk 

Wp ran write a rollprtivp operator for the Ad- 
atom ayatpm aa that for the firat N— 1 atoma plua 


that fur thp Mh at >ni, 

/(:(Ar)./<:(,v- ( 16 ) 

UalnK thp Dirkp atatpa ronatrurtpd by thp repeated 
application of tqa (13), wp ran redurr the matria 
elpment of a rollertive operator of A atoma to 
that of A/-I atoma very eaally. Let ua work out 
the followinK example aa an illuatratiun; 


< r - 1. Ml* I; r; • • • r; - J I /f‘r(A^)l w; r, • • • r,., r - i) 

■ (I 2r){((r ♦»!)(» - Ml* l)|‘^(r- i.Mi* 1 - *';.,1A*’(\ - l)lr - 1, m- i; r,- • • r,.,; 

-|(r-Mi)(r ♦»«* l)|‘^<r- J.Mi* I ♦ i; r', • • • »; .,|«f (.V - l)| - J, m r, • • • » ,.,) 

*|(r»Mi)(r*M;-i)|‘^(r-|,Mi»i:r;-*T;.,|r-l,wt J;r,-**r,.,)p“*^’'‘'‘’^*}. 

( 17 ) 


We can uae the orthonormal property of the Dirke 
atatea to expreaa the laat term In angular brarketa 
In the above equatUm aa the product of Kronerker 

a*a: 

( r - i. Ml* i; »•; • • • r;.,| r - J. mi* i; r, • • • r,.,> 

- •'',)«(»•,. *•;)••' r;.,) 




where 

*(»•,.»•;)• 


I I 
|0 


for r, »r; 
for r, # rj ' 


(IB) 

(19) 


We can alao uae Eqa. (B) to factorize the matrix 
elementa In Eq. (17) and drop common factura on 


both Bidea ; then thia equation will reduce to 
(r - 1; r; • • • z;.,r - i|ffJ(M| r; r, • • • r, .,r - i) 

•d 2*)|(r-l;r;**-r;.,|/fhv-l) 

ThIa la a recuralon formula for one of the RME. 
Wp notice that mduea not appear In thia formula. 

Following Bimllar prncedurea and keeping In 
mind that the aelecllon rulea r, - r', ■ * 1 or 0 
appiv ‘o any of the (r, , r\) palm, we can obtain 
a act c( all the nonvanlahlng recuralon formulaa 
like E |. (20) for HME. To aave apace we will 
drop ihe c, • • • r,., and •'k^ parta of the 

Indlcea. 


(r;r-J|f#;(A0|r:> -i)«|(2» -1) 2» |(r - il/(:(A - 1)| r - i) ♦ (1 2»-)6,.,i-"‘^*“'”^*. (21c) 

(r;r i.l|/<:(.V)|r;r ♦J)-l(2r*3)/(2» *2)|(< ♦ J| - 1)1 r ♦ i) - 1 1 (2r *2)|6,.,e*'''-^i"**. (21b) 

(I : r * i| «:(M| r ♦ i) . - I r(r ♦ l)| -‘V * i| /<:(A- 1)1 r ♦ i). (21c) 

(ri-l:r ♦il«:(V)|r;r fj).( 1 (2i ♦2'l(r l)lr ♦i)-! 1 (2» ♦2)|fl,.,i’’"**'i''^». (21d) 

(r*l;r*i|/f:(;^|i;i--i).|r (r 1)|‘^(» ♦ il A^.V- 1)|»’ - i). (21e) 

(r ♦ 1; ♦ J| A:(.V)1 1 ; ♦ i) - ((I- ♦ 2)/(r ♦ l)]'^(r v || Ar(V - 1)| r ♦ i). (21f) 

(I -l;r-i|K:(,V)|r;r-i).|l 2r |(r - i| 1 )| r - i) - ( 1 2< (2H') 

(r - 1; r - i| A^.V)1 r; r - i) -((r - l)/i |' V - il A|t(.V - 1)1 r - i). (21e') 

(r - 1; r - J|AhM|r; r ♦ i) -|(r ♦ D/rpV ♦ il «i(A- D| c ♦ i). (21f' ) 


( 1 kalualhin ul ( 

From Eqa. (B), we can aee that the evaluation 
of the matrix elementa of collective nperatora 

la reduced to the evaluation of the RME; which. 


in turn, la reduced to the evaluation of the coef* 
ficlenta C , aa can be aeen from Eqa. (9). An 
examination of Eqa. (21) leada to the following 
conclualona: 

(I) The recuralon formulaa can be generalized 
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b>' irpUrIntt .V and r by i and r„ n-apfrtivHy, 

with .*1,2 S, and ran br aummaritMl by 

onr K<^<*ral (ormu.a 

*i/lr r'.- r r' Id 

♦ Ill'll. ' (• ' I .|N, .|r 

( 22 ) 

By rcpfalfd appltcatinn u( thia Krneral formula, 
we ran obtain a Kt'xt'ral rxprraalon for C, aa 


^ ♦'i -I ) 

m •“ » 


(II) The 0|., in Eq. (23) Impilea that the rorf- 
(Iclvnt C!,1f vaniahra unlraa thr a and a' arc 
aurh that whm rxpr<>aa«*d In thr notation of Eq. 

(14) give tf «r) for all / ■ I, 2 up to at Iraat 

;«i-l. In other worda. If we have 


• 0 forf«l.2 f-1 


I #0 fory.f. 


then we will have 


( #0 for I * / and r' • r or r 1 1 
• .0 otherwlar ’ 


Thia intereatInK reault la. of rourae, a ronae- 
qut*nre of the RenraloKlral ronstrurllon. 

(lit) The two funrtiona ! and n in Eq. (22) are 
thr .10 (or txith operatora Hx and H\, ao we 
have 




(26) 


(Iv) II are rrplare (he t In Kqa. (21d')-(2l( ) by 
r *l ami then take lirrmitlan ronjuKatr. they will 
be e.Xartly thr aamr aa Eqa. (2ld)>(2l(|. Thia 
kind of arKument Irada to thr (ollowintt rrlation: 

(27) 

Howrver, uaintt Eq. (23) to rvaluatr C,*7' lor 
larKe ayatrma may atlll be rather tedtoua. In 
See. HID, we will introdure a dUKrammatir trrh* 
nique whirh will help ua to vlaualize thr aituationa 
murh more rlea 'ly. 


II liiaiH'hiii|«li4aiim rrprrwiilalMMi 

The branrhInK dianram* la uaually uard to dr> 
(ermine thr dimrnalon of an Irreduriblo reprr* 
aentation of a permutaliim Kroup % rorreapundintt 
to a binary partition of .V. When thia diaKram la 
uaed In the dlaruaalon of a ayatem of eirrtrona, 
the abariaaa and the ordinate rrprraent thr num- 
lor of eirrtrona and thr total apln. rraprrtively. 

In our raae, they will reprearnt thr number of 
atoma and thr rooperatlon number, rraprrtively. 
We will represent a Dirke state |r, »s;r, • • • r,.,; 



T«oN 

> * p I ' < n I 


ITU. 1. Graphical repreaentatlon of an example of reduced matrix element ( jt llilili2| l*i|ltU of a collective 

operator of a nine-atom ayatem corresponding to an arbitrary mode k. The Iwckgrouixl Is the branching diagram usually 
used to determine the dimensions of Irreducible representations of permutation groups S/, corresponding to binary parti- 
tions of .V. The abscissa denotes the numlier of atoma In the system or a subsystem, and the ordinate represents the 
cooperation numlier. The letters at the bottom Indicate the types of sections, defined In Fig. 2, of this example d las ram 
which Is of type .VC3. 
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by a alK'EaK Hnr un thia branrhlna diaKrant, and 
a pair <j( aurh linra will rrprrarnt a rvdurrd tna* 
IrU rirmrni aa ahuwn In Fig. 1. In Uila diagram 
rrprramUtlon, Ihrrr la nu way lu tell the m valur 
uf an rigcnatalf. Fortunalrly, a rMiured matrix 
rirmmi la lndrp<-ndrni ol m; au wr npt<d not worry 
about It. Brrauar o( Eqa. (26) and (27). wr alao 
nrrd nut dtatlnKulali wtirthrr tlir operator 1a M*C 
or K] and witirli line repreaenla tlie atate im the 
right lA the opera*'"’ the ket vector, and which 
lln ’’epreaee' /•late on the left, the bra vec* 
tor. 

From Eqa. (21) and (27), we can aee that the 
reduction of a RME uf N aluma to that uf S- I 
atoma can occur In alx different alluallona. They 
are repreaented by a aecllon ol llnea, which rep- 
reaenla thta PME. between two vertical llnea 
curreapunding to .V and .V - I, The alx different 
aeclluna are ahown In Fig. 2. The factora 
Mr„r;;r,.„r;.,) and /f(r„ r',; r,.„ r',.,) can alao 



KIG, 2. Six dlitinct type! of aectlona o; the graphical 
repreaentatlon of a rrducnt matrix element of a collec- 
tive operator. 


be ctNialdered aa funcllona uf tne type at a aec- 
tlun and the r value uf one of the four cumera, 
aay the upper-left one, of Ihia aeclltm auch aa 
/(/t,r) and g(/f, r). Theae funcinxia are Haled In 
Table I. 

I I «mifl 

With Fig. 2 and Table I aa gulden, we ran now 
eaally uae Eq. (23) to evaluale, aa an llluatratlon, 
all the ruefficlenta of the RME of nine atoma 
repreaented by the diagram In Fig. 1 aa folluwa: 
The firai nunvanlahing value of (r, -r') uccura 
at I ■ 4; ao we have 

C,*C, 

witere, aa well aa In the following, we liave dropi'ed 
the other aubncripta and auperacripta In the nu- 
tation of the coefficienta. 

A rcmmon factor of all the r>'mvanlahtng coef- 
ficient! la 

/ 4-, • y (*•. I )/(f. I)/ (A. I ) I (/>. 1)/ (F. 2) 

-(I /J)|-2 (4«2)‘*|(|Hi)^f • -1^ /f-. 
TTien we have 

t‘4 ’HU*. 4^ •( - i)(-l^)*T' 

f, - v(/f. 1 ).f i )/ ,^ . ( - iMi M - i^)*f . 
f, -ifU. i)/ («. D fiH. i) f 4 - • (i )(f Mi )( - . 
C , ./f(A. 0) f (A. i )/ («. I )/ {U, i )y 

So the (Inal reault la 


T\. l.r I. Parlors In tite recursion formulas for dif- 
ferent tvp4’S of scrtlona. us defined In pig. 2 
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IV. VSVMriOIM RIM IIS MM I \M«.t V 

From the phyairal point ul view, we are mure 
Intereated In the aaymptutlr reaultn fur Urxe N 
(* 10*”) than In the exart reaulta for amall .V, When 
A' la ao lante, we have r»l for an uverwhelmliiK 
ma)orlly of raarw; no we can aaaume r to be of 
the name order an A without intr<xfuclnK any aer- 
loua errorn. 

From Table I. we nee that In the aaymptotic 
limit tne / factor la I for a nectiun of type A, U, h, 
or F and of the order l/A' for a aectlon of type 
I or O, Therefore, In the diattram repreaenlinK 
a RMK, whenever a aectlon of type C or U uccura, 
the value of thin RME will be reduced by a factor 
of the «irder of I A’. For the convenlt*nce of atate* 
menta, we nhall refer to a aectbai of type /> or C 
an branchlnic or rejolnlnic aectlon. reapecllvely. 

K I .ipiiliifi ji t lawirHiinon irf irducml m*int rirmeiilt 


it{S~d - II. Therefore, the condition that we 
can diacard proceaaea of type C2 In U(A-J - 1) 
which implien 

d/W-1. 

On the other hand, for a fixed ket vector 
(r, r,.,), there arc «f poaaible bra vec- 
tora<r ♦ 1. »n» 1; rj* • • ri.J and W-<f - I p<MBible 
bra vectora - 1 .»m» l;r;** •»;.,! to form pro- 
ceaaea of type .VL 1. So under condition (29). we 
can diacard proceaaea In which r Increaaea. 

Cimaider a diagram of type .VCl with the branch- 
init p«)int uccurrlnx at the A,th aectlon auch that 
1 A'l - »d,. There will Im- «f,(W, -if, - I) poaaible 

waya to form the “liHip" of a diaitram of type At 3 
correapondlmt to thia fixed diaitram of type At I. 
Therefore, the condition that we can diacard pro- 
ceaaea of type At 3 la d,(.V, -if, - l)«* ,V{, which 
impilea 

d, A,«l. (29- I 


We can claaalfy the RME accordintt to the topo- 
loitlcal typc‘B of their diaRrama, an ahown In Table 
II. Obvl II ily, type t'O, repreaentlnii a diaitonal 
element, la the moat Important one amonR the 
r-connervInR prociwaea; and the name la true 
with type .VC' I amuni; the r-nonconservlnit proceaa- 
ea. The value of each of an RME of type C2 or 
At'3 la amaller than the leadlnR one by a factor 
of the order of I /A*. Therefore. If the number of 
poaatble diaitrama of theae typea correapondlmt 
to each of the leadlnit one la mi'Ch lean than .V*. 
then we will be aafe to diacard them. 

Let UB define 

if a i .V - r, (28) 

Then for an eiitenalate |r, we will 

have r,,, - r, • - J exactly d timea. In other 
worda, the diagram repreaentlmt thia atate will 
drop exactly d timea. 

Fur a fixed ket vector I r, mi; r, • • • r*.,), there 
la only one bra vector to form a proceaa of type 
CO; and there are lean than if(.V -if - 1) poaatble 
bra vectora to form proceaaea of type C2. The 
btter can be aeen aa followa: In the branching 
diagram, one can aee eaally that one end of the 
"loop" In a diaRram of type C2 muat occur at one 
of the d aectlona with r ,, , - r, ■ - J, while the 
other end can occur at any of the other A’ -if - I 
aectlona with r,., -r, ■ So there are 
(f(.V-if - 1) poaalbtittlea. However, aome of theae 
poaalbllttlea may have to be dropped becauae of 
the reatrli Ton that r, can never be nexattve. 
Hence, the real Doaalbtlltlea muat be leaa than 


Fur an overwhelmimc majority of caaea, we can 
nay that condition (29) impIlea rondilion (29 ). 

Fortunately, from the jihyBtcal point of view, 
we are moat Intereated In altuatlonn where con- 
dition (29) la aatlafled. 


TARI.E II. (’Uaalflrallon of reduced matrix elementa 
aiTordlna to ilip to|M>lo|[lcal tvpea of thrtr diaarama and 
the order of maanitude of each type. 
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Am ran hr nrrn (rom Tahir II, a dlaxram al 
KME ronalMla o( MinRie Unrn and douhlr tln«’a 
altrrnaMvrly ronnrrtrd hy hranrhlnR and rrjoInlnR 
arctlona. Wr Mhall rrirr to a alnxlr llnr aa an 
$ "propaRaior" and a douhlr llm* aa a T "propa* 
Ralor.” Howrvrr, the flral almtlr llnr mual hr 
Irratrd dKfrrrnlly and ahall hr rallrd Ihr "Km* 
rraior’* lx‘cauar, while othrra rontrihutr only a 
common farlor to all the nonvanlahlnR corKIrlrnla 
Ihla alnRlr llnr dricrmim'a the alrurturr 
or variation of the roefficirnta. We ahall denote 
the contribution to the RME by an ^ propaxalor, 
a T propaRator, and the Renerator hy ii, T, and 
V, rraprctlvrly. 

/ .V frnfugtl'in 

An ^ propuRator la a conaecutive product of a 
mixture of /(A, r) and /{B, r) factora. It turna 
out that the reault of thla product alwaya Rlvea 
the numerator of the fir at factor divided by the 
denominator of the laat facto” Thla can hr aeen 
eaaily by conaldeiinR the * ' • ir aimple ex- 

ample: 


y <A, r, >/(A, r,. ,)/ (ff, (A. r,.,l 

, 2r, ♦ I 2r, ♦ 2 2», ♦ I 2r, ♦ 2 

2r, ♦ I 2r,* 2 2r, ♦ I 2r, 4 2 2r, * 3 

2n ^ 2r. 

’ 2r, ♦ 3 " 2r,., ♦ 1 

Therefore, an ^ propaRator la a function of the 
typ«>a and r valuea of the hrRlnnlnx and endinR 
ai-ctlona only. And wr have 

i(A,».;A,r,)"i(A,r,; W, r^)"2r, (2r^ ♦ lit 

(30) 

i(ff, r,;A, rf)"S{B, r,; B, (2r, ♦ 2) (2r, ♦ 1) . 

For larRe valuea of r„ theae two expreaalona are 
almoal the aame and can be aummarlzrd aa 

r, , (31) 

2 1 fntptfmlim 

A T propaRator la a conaecutive product of a 
mixture of J (£, r) and t (/*, r ) factora. Let ua 
firat obtain the reault of a conaecutive product 
of pure ./ (£, r) factora. We have 


/ (A‘. ZV )/ i‘- , I ,)/{£■, r, ,,)•• •; (£, »> ., )y ( A, r,. , ); (A, r , ) 


Then conaldrr the case when one al the factors 
J (A, Kf ) Is replaced by / (A, r, ). We have 

/ (A, r, )• • •/(£, r, .,)/(£', r, \f (A, r, , ,)• • •/ (A, r,) 

_ ( {2r,- l)2r, 2r, ♦ 1 (2r, -2><?r, -1) V'» 

“\(2r,-l)2r, 2r, - 1 2r,(2r,*.' ) 

-li(r,)|(2r,- l)(2r,) (2r,)(2r, 4 1 ))•'•, (33) 

where 

*(r, ) a 1 (2r, - 2 )(2r, 4 1 )/(2*j - 1 )(2r, )) *'* 

-l-l2r,)'’4-...l-o(iV-'). (34) 

We see that the deviation of Eq. (33) (rom Eq. (32) 
Is of the order of .V*. When there are il' t (A, r) 
(actors, the deviation will be of the order of 
d' .V’. Therefore, for larRe values of r,, we have 

r, . (35) 

It is InterestlnR to note that 5 and T propaRators 
ire the same asymptotically. 

( liijsiinal mlumi matnt rkniirnlt and ihr frorralnr 

A dlaRonal HME is represented by a diaRram 
of type CO; i.e., one sinRle line only. Suppose 


1 2r. 2r. 4 1 2r,- 3 ^ 2»,- 2 . 2»,- 1 

, 4 1 2V| 4 2 2r, 4 3 2r^- 1 2r, 2r, * 1 / 

lH2r,) (2r,H2rj* 1)1’'’ . (32) 

that sections of type A occur S-d - 1 times and 
those of type B occur (/times at the (i,- l)th 
section with a ■ 1, 2, ... , d. AssuminR that the 
last section is of type A, we can calculate the 
coefficients (or (our different cases accord- 
InR to the types of the (/- llth and the ith sections, 
namely yt/l. A//, BA, or AA, as follows: 

CJ*^ ■ A , r, . , )S( A , r,; A , r - i ) 

■(1 2r,)(2r« 2r)-l 2r, ,3(>a) 

*'(,.i)Ji(A, r,,.,;A, r -i) 

-ll/(2r,^4 l)Jl(2r,^4 3) 2r) 

• (1 2 r ) l ( 2»-,^4 3 )/( 2 r ,,4 D ) , ( 36 b ) 

C ■ >{(A, r,,.,)5(A, fi. ;A, - J) 

-I- l/(2r,.4 2))|2r,^ 2r) 

■-(1 2» )l(2r,,)/(2r,, 4 2)1, 

Cf/ - A, r,,. , )S( A, r,. ; A , r - i ) 

■[-1 (2r,, 4 2)ll(2»-,, 4 2) 2r) 

■-1 2r, (36d) 

when* we have dropped th< reRular subscripts and 


» 
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«uprrs''np(ii of Uie ro«fflrlpnU. 

The probahURv that AA will occur Is about W that of AB or BA la about d(N -<fl.V: and that of 

BB la about (f'/AI*. Therefore, under condition (29). casea of AA are an overwhelmltiff majority and thoae 
of BB are neKlUi^'*?- 

Let ua conalder a conaecutive pair of AB and BA, Thla pair will contribute the followlnft to the HME: 

^ Cf* e " •'* 

■(1 ♦ U r,, 1(1 ♦ 


The value of P, aa Kiven by Eq. (37) dependa on 
the Index (K -lc,)*(X,^-X,,.,); which, In turn, 
dependa very critically on how we order the atoma 
In the aample. Aaaumlng that the atoma are 
randomly dlatrlbuted, we can order them in auch 
a way that two conaecutive atoma are very cloae 
nel|(hbora and (X,^-j{|,.,) la almoat perpendicular 
to (i( -!(,), aa illuatrated in Fig. 3. Ihippoae In 
‘Mb way we can manage to reduce the Index to a 
value of the order 

(K-U,)«(5?,,-X,,,,)-0(t); (38) 

then we have 

P,mU(€ N)*U(1 N), (39j 

where the laat atep la taken becauae we believe 
that the ordet of c la liifther than that of I S. 

Now we can write 

(r;al/f< |r;a)» 


where a poaaible error of the order of d^/S' la 
Introduced becauae we have Ignored the exiatence 
of BB, The prime on the summation slRn indicates 
that i, - 1 and i, are to be skipped; this means 
that the number of terms to be summed is exactly 
2r. This would be true even If we had Uxen cases 
of BB into consideration because for every term 
of BB, a term of AA is canceled almost exactly. 

If the last S' ction of the diaKram Is of type B, 
then the last factors in Eqs. (36) should be 
/{B, r instead of J (A, r -i), and 2r In all ex- 
pressions In this subsection should be replaced 
by 2r * 2. But this is of no significance U all. 

The generator of an RME Is represented by the 
first single line of its diagram. Suppose this 
single line consists of A', - 1 sections with 
A', - <f, - 1 of them l>elng cf type A and d, of type B. 


the diagonal RME as 

M 


2r 2- ^ 
I • I 


■ 

♦ ^ P, + 0(d‘ iV’) 


The expression for the generator can then be ob- 
tained from ttiat for the diagonal RME In Eq. (40) 
by replacing S and r by A, and r,,^, respectively, 
as follows: 

i> 

57 — ♦D(d,t A,). (41) 

!•! 


II On-dumunol rrdiH'«d nulrix rImirnU 

With the expressions for propagators and the 
genera'.or available, we can very easily wrlie 
down the RME for all types of diagrams. For 
off -diagonal RME, the generator la always fol- 
io v.d by a section of type li. It Is convenient to 
combine the contribution of these two to the RME 
once for aU ac follows: 



KIG. 3. An illuRtration of the way to choose the se- 
quence of atoms in order to reduce the possible error 
introduced In simplifying the expression for a diagonal 
reduced matrix element or a “generator.” 
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<n) Tyfif NCI. 

(r*l;a'|«rk;3)-«^‘(.V,)T(r,,.„r,.,). 5ir(^ ^ , _ ,43, 

*1 I • I 

(hi Type C(MI. A dlaKram of type C'(2i>) consUta of the K«'nerator, n T propaKatora, and n S propaxatora 
connected by a aectlona of type U and n aectiona of type C. We can write the RML' aa 

(r;o'|/fi|r;0t)-vfLV,)7(r», )/(G, r,, r,,)* • ‘Mr,,,,,, r,.,) 

fr> Type NC(Jh ♦ ll, A dlaxram of type A^C (2n + 1) conalata of the x^nerator, (a ♦ I) T propaxatora, and 
a ^ propaxatora connected by (a-f I) aecttona of type U and a aectlona of type i'. We can write the HME aa 

2Kifii:^ 5 v (si: t ■ |«| 


V. t Vl INDRK AL S\Mm AM) 1 RANSII K)N 
PRUBARILIIIfS 

To obtain numerical values of HME by ualnx 
Eqa. (40), (43), (44), and (43), we replace the 
summation by Intexratlon. The result of such 
Intexration will, of course, depend on the boundary 
of the sample. This Is where the sample-shapi' 
effect comes in. The most popular sample shape 
used In the discussion of superradiance has been 
the circular cylinder.* In the literature, cases 
of spherical'* and rectanxular sha|X'" have also 
been considered. In this paper, we will cor.slder 
the case of a circular cylinder only. 

A ( >lindrical vaniplr 

Consider a sample of atoms confined in a space 
of the shape of a circular cylinder of length / and 
radius a. We can express the summation in Eq. 
(43) as an integration over the sample space as 
follows; 

I ' I 

■ —7- f </« 

MllU Ja 




where (k -k,)<,p are the transversal and axial 
components of (k -il,), and the a axis has been 
chosen in the k, direction. Carrying out the Inte- 


gration in Eq. (46), we obtain 


slnl(it-it.)W 2j 2J,((it-ic,)pu) 


l(R-k,),/ 2 ] 


(47) 


where J, is the Bessel function of order 1. 
Now we can write Eq. (43) as 

ir± l,a'|«t|r,o) 

-(1 2rKf''.(E. 


where we have assumed that the order of the atomt 
is such that a disk of radius a and thickness dz 
is completely filled before we go to the next disk 
along the axis of the cylinder; we have also de- 
fined 

/,s/(,V, N). (49) 


B T'inUliiin probahililirs 

With our knowledge of the matrix elements but 
without solving' the master equation, we should bt‘ 
a)<le to get some Idea about the time evolution of 
tie atomic s'stem. We assume that at time (■O, 
tl.7 system consists of completely inverted atoms 
w.'th r ■ j.V and a vacuum of the radiation field. 
The initial evolution when the value of r is still 


i i \ i i 


12 


DIAGHAMMATIC TECHNIQUE FOR CALCULATING MATRIX 


585 


cloM to Will b« esaentlally determined by 
the collective oper?* r Rf . Suppose at a certain 
Instant the atomic system Is represented by an 
eigenstate |r, m;a); we will consider the transl* 


tlon p'obabllltles to the slnRle state |r, m- l;a) 
and to the whole set of stales |r - I, m - 1; o') , 
respectively. They will be proportional to the 
(ollowlnx quantities; 


(a) r-i'onuerriHii prucfMu. 

l;a|Rf |r,m;a)|*«(r l)|(r; a |Rf |r; a)|* 

• (r ♦ m)ir -m* 1) 

(ht r~mmctmserviHii prucemtes. 

« 

^ |<r- 1, w- l;a*\Hf\r,m;a)\*‘(r*m){r*m - l)y2'l<»'- 1; k; «)|* 


sin 

T 


(S-K,),/ 2|* • 


(50) 


* (r ♦ m)(r ♦ m - 1 ) di, f „ a) - « ' I * '(2r F 

"t** 

(r ♦ m)(r ♦ >M- 1 ) / , sln''l(k - jl.). / 2| » .... 

Tr 1‘- [(r-R,).r2l^ [(d-M' ‘ * 


From Eqs. (50) and (51) we can see the general 
behavior of the evolution of tltc atomic system. 

At the befUnnlnfi, when m la close to r, the fac- 
tors (r + »i)(r-»w* l)of Eq. (50) and (y+w)(r-*-»M-l) 
2r of Eq. (51) are of the same order of magni- 
tude. So the relative Importance of the two lUnds 
of processes Is determined by the other factors; 
namely, |F(ic, l,a)\* of Eq. (50) and 1 -|F(ic, /,«)|* 
of Eq. (51). Rehler and Eberly* have calculated 
numerically the function |F(k, f,a)|', ''^..ch shows 
a high peak along the direction of i(,. Therefore, 
for radiation around the direction of k„ It is 
predominantly r conserving; and for radiation 
in other di>ections, It is predominantly r non- 
conserving. Later on, as m approaches zero, 
we liave (r * m)(r -m* 1)» (r + m)(r-*- wt-1) 2r, 

If r» 1. Then the r-conservlng process dom- 
inates, and the radiation will be confined within 
a small solid aiiglc around the direction of H,. 

We can see then the condition for superradiance 
is that the st te of low m must be reached before 
r orops too low due to r -nonconserving processes. 

VI SUMMARY 

In the Dicke-Schwendimann approach to the 
problem of N two-level atoms spread over large 


regions, interacting with multimode-radiation 
field, the ,na)or difficulty has been lack of ex- 
plicit expression for the reduced matrix elements 
tRME’s) of collective operators which determine 
the static aspect of the problem. The aim of this 
paper has been to clear up tliis difficulty. The 
genealogical construction was adopted to obtain 
Dicke states corresponding to a specified mode. 
Then a set of recursion formulas of the RME’s of 
collective operators corrcspo.idlng to arbitrary 
modes was obtained. Using the branching diagram 
usually used to determine the dimensions of bi- 
nary Irreducible representations if the permuta- 
tion group S„, a simple and systematic method 
was developed for calculating RME’s for any 
number of atoms in a system. However, from 
the physical point of view, one is most interested 
in the asymptotic resultsfor large Af. In this asymp- 
totic limit, it Is found that the order of iiiagnitude 
of an RME Is determined by the topological struc- 
ture of Its diagram. Explicit expressions have 
been obtained for RME's with diagrams of ail types 
of topological structures. They are presented In 
Eqs. (40), (43), (44), and (45). Equation (40) gives 
the diagonal r-conserving RME with r close to 
which is not much different, especially wiien the 
summation Is replaced by an integration, from 
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the case with r ■ J.V which Is well known. Mow* 
ever, because this expression has been obtained 
with the order Ol maanitude of possible error 
clearly In mind, the result certainly enhances our 
confidence In extendlna the usaae of this expres* 
slon from the case of r ■ ^ .V to cases of r close to 
i.V. Equation (43) Rives the most Important off* 
dlkRonal and r-nonconscrvlnp ItME. This Is, per* 
haps, a new result, althouRh .Schwendlmann has 


speculated about a similar expression (Ref. 4, p. 
276, footnote). Anvway, In my Judnement, the 
most Important contribution c»f this (taper Is In 
KalnlnR InsiRht Into r-nonconservlnR processes. 
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